Nickel assisted healing of defective graphene 
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The healing of graphene grown from a metallic substrate is investigated using tight-binding Monte 
Carlo simulations. At temperatures (ranging from 1000 to 2500 K), an isolated graphene sheet 
can anneal a large number of defects suggesting that their healing are thermally activated. We 
show that in presence of a nickel substrate we obtain a perfect graphene layer. The nickel-carbon 
chemical bonds keep breaking and reforming around defected carbon zones, providing a direct 
interaction, necessary for the healing. Thus, the action of Ni atoms is found to play a key role 
in the reconstruction of the graphene sheet by annealing defects. 

PACS numbers: 61.48.Gh,61.72.Bb,68.65.Pq,81.16.Hc 



I. INTRODUCTION 

Carbon structures — graphene and nanotubes in par- 
ticular — are the focus of substantial research activity, 
due to their promising electrical, mechanical and opti- 
cal properties^ Particular attention has been paid re- 
cently to single layers of graphite, i.e. to graphene, 
and to its novel electronic properties^ In many appli- 
cations, optimal performance requires the control of its 
structural properties, which remains a significant diffi- 
culty. It is known for a long time that Catalytic Chemi- 
cal Vapor Deposition (CCVD) of hydrocarbons on reac- 
tive nickel or transition-metal-carbide surfaces can pro- 
duce thin graphitic layers. 5 However, for the most part, 
the large amounts of carbon atoms absorbed on nickel 
foils formed thick graphite crystals instead of graphene 
films. Few layer graphene is now commonly produced by 
CCVD at temperatures ranging between 350 and 1000° C, 
using metallic substrates (such as Co, Ni, Ir, Ru) that 
catalyze the decomposition of the carbon-bearing gas 
precursors and make the growth of carbon nanostruc- 
tures^ — Large- area graphene films of the order of cen- 
timeters are now grown on copper substrates possible*^ 

Despite this considerable progress in synthesis pro- 
cedures, the details of the microscopic mechanisms in- 
volved in the growth of carbon structures are still lack- 
ing. Atomic scale investigations are far from being fea- 
sible experimentally, whereas computer simulations al- 
low such investigations. We have recently developed a 
tight-binding (TB) model for nickel and carbon that uses 
Monte Carlo simulations in the grand canonical ensem- 
ble (GCMC) to study the formation of graphene from 
a metallic substrate j 14 i 15 The GCMC method works in 
an open system and can simulate the growth of a car- 
bon structure under conditions of fixed temperature and 
carbon chemical potential, closely reproducing the exper- 
imental conditions where the carbon chemical potential 
is fixed by the catalytic decomposition of the carbon- 
bearing gas. Hence, we were able to understand and 
explain the nucleation of graphene on Ni (111). The Ni 
surface catalyses the nucleation process: carbon atoms 



are confined in its surroundings until a critical surface 
concentration is reached. Carbon chains appear first, in- 
teracting weakly with the surface and becoming more 
and more stable; then nucleation of graphene takes place. 
The energetics of the carbon-metal interaction, a lim- 
ited carbon solubility, and a tendency for expelling the 
sp 2 carbon structures formed at the catalyst surface are 
therefore important. 

However, for a large range of chemical potentials, 
we observe the presence of defects in the grown sp 2 
structure. Likewise, in the case of carbon nanotubes 
grown from a metallic nanoparticle, whether the em- 
ployed method is empirical, 16 or semi-empirical j 17 i 18 all 
the final configurations are plagued by a high concen- 
tration of atomic-scale defects. These include, but are 
not limited to, heptagon-pentagon topological defects, 
adatoms, and atomic vacancies. Furthermore, the esti- 
mated experimental growth rates are in the nanometer to 
micrometer per second range j 19 i 2Q Such scales are largely 
beyond the capabilities of numerical simulations, mean- 
ing that only some elementary steps can be studied, or 
that the growth conditions imposed in the simulations 
are orders of magnitude too fast, leading to very defec- 
tive sp 2 carbon structures. In this context, several groups 
have sought to characterize the crystallization of carbon 
network using numerical simulations ] 17 i 21 i 22 

In the present work, we investigate the healing pro- 
cesses of defective carbon structures at atomic scale. We 
have developed an improved Tight-Binding Monte Carlo 
(TBMC) code that allows us to perform systematic and 
long Monte Carlo runs. We use it to study the evolution 
at finite temperature of graphitic structures previously 
grown by grand canonical Monte Carlo (GCMC) simula- 
tions. We also investigate the role played by the catalyst 
metal in the healing process of graphitic structures. 



II. METHOD 

We have developed and carefully tested a model for Ni 
and C, in a tight-binding framework j 14 i 15 The total en- 
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ergy is a sum of local terms: an empirical repulsive term 
and a band structure contribution including s and p elec- 
trons of C, and d electrons of Ni. Local densities of elec- 
tronic states are calculated using the recursion method. 
Only the first four continued fraction coefficients (a n , b n ) 
are calculated exactly. The continued fraction is then ex- 
panded to the N th level using constant coefficients equal 
to d2 and 62 (here N = 40). This approach provides us 
already with a good description of the angular contribu- 
tions to the energy and leads to a relatively fast scheme. 
This atomic interaction model is then implemented in a 
Monte Carlo code, based on the Metropolis algorithm, 
which allows for structure relaxation. To go further, we 
have improved the Monte Carlo code efficiency in canon- 
ical ensemble by modifying the algorithms. 24 Switching 
from the standard recursion algorithm to a direct calcula- 
tion of the moments is possible and more efficient when 
only four moments are considered. Indeed in a Monte 
Carlo algorithm, we can calculate only the moments con- 
tribution that have changed at each step. A second point 
concerns the way to calculate the local densities of states, 
by using an analytic integration of the continued frac- 
tion expansion^ This option is faster but is less stable 
and more difficult to implement. By combining both de- 
velopments, we have obtained a faster code, that makes 
the exploration of longer time scales and more extensive 
investigations possible. More details will be given else- 
where. We have studied the influence of temperature 
variations on defective graphene by using this approach. 
The idea is to start from graphene sheets containing de- 
fects and submit them to high temperatures up to 2500 
K, which is a typical temperature used in graphitization 
processes^ The system is then able to overcome high 
energy barriers and reach new equilibrium states corre- 
sponding to less defective carbon structures. 



III. STONE- WALES DEFECT 

As a first step towards a complete understanding of 
the healing process in graphene, we study the influence 
of temperature variations on the Stone- Wales defect. The 
Stone- Wales (SW) defect is a 90° rotation of a C-C bond 
in the hexagonal network with respect to the midpoint of 
the bond. This leads to the formation of two pentagons 
and two heptagons, replacing four hexagons. This de- 
fect, studied extensively, is believed to play an important 
role in the structure transformation of different carbon 
nanostructures, such as the coalescence of fullerenes or 
the melting of carbon nanotubes i 27 i 28 According to ab 
initio studies, its energy is of the order of 5-6 eV with 
an energy barrier energy of the order of 10 eV^ — 

Without the presence of adatoms which can reduce 
this value within the range 0.7—2.3 eV^ 2 - it is then 
impossible in practice to observe the formation of SW 
defects at equilibrium. The barrier for healing of an 
existing SW defect on the other hand is obviously much 
lower, of the order of 4 eV or less. To go further, we 
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FIG. 1: (a) Reaction pathway for Stone- Wales transformation 
observed during a MC simulation at 1000 K. (b) Schematic 
qualitative representation of the activation barrier for healing 
the Stone- Wales at high temperatures. 



use Monte Carlo simulations in the canonical ensem- 
ble to investigate the behavior of this defect at high 
temperatures and without constraining the system. 
Indeed, many ab initio calculations up to now have been 
performed at K by imposing an in plane Stone Wales 
— graphene transformation. In the present work, the 
system consists of a central defect site within a supercell 
sufficiently large to minimize boundary effects on the 
energies of interest. We adopt a structure containing 
500 atoms. The full structure is relaxed via standard 
MC displacements steps. FigJIJa) shows that when 
the defect is subjected to high temperatures (> 1000 
K) the median bond separating the pentagons and the 
heptagons completes a 90° rotation which results in a 
perfect sp 2 hexagonal geometry. During the simulation, 
we have observed that this thermally activated healing 
process involves an out-of plane path (FigJIJb)) (see 
also Ref . (3lj ) . This reaction pathway seems therefore to 
reduce significantly the activation barrier. 



IV. DEFECTIVE GRAPHENE 

We now consider a more complicated situation. A 
graphene sheet grown by previous GCMC simulations^ 
is submitted to temperatures ranging from 1000 to 2500 
K. The initial configuration contains 72 carbon atoms 
where 70% of the sheet is populated with defects (see Ta- 
ble H]) . The convergence of the total energy as a function 
of Monte Carlo steps is controlled and compared with 
variation of the number of rings in the system. The sim- 
ulation is stopped when the total energy no longer varies 
on the average, which implies that the system has reached 
a Gibbs energy minimum. Fig |2] shows the final states of 
the MC runs for four temperatures (T = 1000, 1500, 2000 
and 2500 K). Table U presents the corresponding types of 
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TABLE I: Number of rings in initial and final configurations 
for different temperatures. 
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FIG. 3: Comparison of the total internal energy with the vari- 
ation of the number of polygons as function of Monte Carlo 
steps in an isolated graphene sheet at (a) 1000 K and (b) 2500 
K. 



FIG. 2: Equilibrium configurations at (a) 1000 K, (b) 1500 
K, (c) 2000 K and (d) 2500 K. 



polygons. At 1000 K, a significant amount of defects re- 
mains in the final configuration. Upon heating, octagons 
and nonagons are the first to be corrected, implying that 
their activation barrier is the lowest. Pentagons and hep- 
tagons are the most difficult to heal and remain in large 
concentration until 2500 K. Consequently, topological de- 
fects of different sizes are thermally activated at different 
temperatures. During the simulation at 1000 K, we ob- 
serve a stepwise decrease of the total internal energy (see 
Fig|3ja)). The steps coincide with a significant increase 
in the number of hexagons. Indeed, the configuration 
that is favored by the carbon atoms is a three-fold en- 
vironment, where the bond lengths and angles of an sp 2 
configuration are respected as in the hexagon. Such a 
geometry greatly stabilizes the system. This tendency is 
further confirmed by calculations done at 1500, 2000 and 
2500 K as shown in Figj3^b). However, the steps are less 
visible at high temperatures (> 1000K) due to the large 
fluctuations of the system which results into noise in the 
plots. As expected, the presence of hexagons greatly sta- 
bilizes the system and tends to heal all other types of 
defects. As seen in FigHJ where the evolution of the sys- 
tem at 2500 K is presented, we observe a predominance 
of hexagons at the end of the simulation. The final con- 
figuration, corresponding to a carbon structure with 29 
hexagons, is indeed an almost perfect graphene sheet. 




FIG. 4: Evolution of the defective graphene sheet at 2500 K. 
The transformation of several defects in an isolated graphene 
sheet, (a) Linear chain of C helps the transformation of var- 
ious defects such as vacancies and large polygons, (b) Va- 
cancy defect, (c)-(d) Stone- Wales defect transformation dur- 
ing heating, (e) Final configuration (top and side view) 



V. ROLE OF NICKEL 

We further investigate the correction and migration of 
defects in sp 2 structures by adding to the previous system 
the Ni(lll) lattice from which the graphene has nucle- 
ated (FigE^a) and (b)). We therefore seek to understand 
the role played by the metal in the healing of defective 
carbon structures. The NiC system, containing 216 Ni 
atoms, was subjected to high temperatures ranging from 
1000 to 2500 K. As in the case of the isolated graphene 
sheet, we notice that at 1000 K, the sheet is far from be- 
ing healed and a non negligible concentration of defects 
remains (FigJS^c)). Once again, large rings are healed at 
1000 and 1500 K whereas pentagons and heptagons can- 
cel out at higher temperatures. At 2500 K the graphene 
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FIG. 5: Graphene sheet in the presence of the Ni lattice, (a) 
Side and top views of the initial configurations. Equilibrium 
configurations at (c) 1000 K and (d) 2500 K. 



sheet is completely healed as seen in Figj5jd). This sug- 
gests that the healing of defects is thermally activated as 
observed for the isolated graphene sheet. By compari- 
son with the case of graphene sheet without catalyst, we 
deduce that Ni plays a crucial role in defect migration. 
First, the Ni lattice has contributed to a more rapid stabi- 
lization of the system: the graphene sheet was corrected 
in less Monte Carlo steps. At 1000 K, for an isolated 
graphene sheet, 2000 MC steps per atom were necessary 
to form 14 hexagonal rings whereas in presence of nickel 
atoms the same healing is obtained with only 600 MC 
steps per atom. Then, the pentagons and heptagons, 
which were extremely difficult to heal at low tempera- 
tures become more reactive in the presence of Ni. To 
understand the role played by Ni atoms in the healing of 
defects, we have investigated the step-like increase during 
the variation of the number of polygons in a simulation 
at 2500 K (FigJBJ). This detailed analysis suggests that 
Ni atoms can interact with the defects of the graphene 
sheet. During the simulation, a defective zone is identi- 
fied: the red one in Figj6j The Ni atoms are seen inter- 
acting precisely with the defective parts of the graphene 
sheet whereas no Ni atoms are found in the proximity of 
the areas of the graphene sheet where a hexagonal geom- 
etry has been established. This may not be so surprising; 
defects change essentially not only the electronic proper- 
ties but also the chemical properties of carbon structures, 
being centers of their chemical activities^ several ab ini- 
tio calculations report on the strong interaction of transi- 
tion metal impurities with defects^ We further observe 
during the simulation that the interacting Ni atoms cat- 
alyze certain reconstructions. The step-like increase in 
the number of polygons corresponds to Ni atoms inter- 



acting with defective areas which help to establish about 
six hexagonal rings (Fig|6jb)). Once again, these obser- 
vations are in good agreement with ab initio calculations 
demonstrating that Stone- Wales transformation occur- 
ring during growth are made easier by the presence of 
additional metal atoms, reducing the barrier from about 
6 eV to 2.80—3.5 eV j 35 i 36 Moreover, recent in situ X- 
ray diffraction measurements of amorphous carbon on Ni 
surfaces have suggested graphitization by direct metal- 
induced crystallisation, confirming the key role played 
by the metallic substrate to heal defected carbon struc- 
tures. 39 
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FIG. 6: Variation of the number of polygons as a function of 
the MC steps at 2500 K . Ni helping the correction of defects 
corresponding to a step-like stabilization of the system. We 
observe that the defective zone (in red) vanishes. This is 
mainly due to the Ni atoms which interact with this area as 
seen on the snapshot. 



VI. CONCLUSION 

In conclusion, based on tight-binding Monte Carlo sim- 
ulations, we have studied the healing processes of carbon 
structures. High temperatures are capable of healing the 
Stone- Wales defect by adopting an out-of-plane path that 
significantly reduces its activation barrier. A study of an 
isolated graphene sheet has demonstrated that defects 
were thermally activated and that high temperatures are 
able to heal the structure. Once this sheet is put in con- 
tact with the Ni lattice, the defects are healed at a faster 
rate and at lower temperatures. We have also observed, 
during our simulations, that the Ni particles participate 
actively in the healing process. Therefore the presence of 
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Ni favors not only the nucleation of the graphitic struc- 
tures, but also the healing of defects. 

Finally, in a more general context the results reported 
here are an illustration that defects in sp 2 carbon struc- 
tures are highly sensitive to the temperature and to the 
presence of a transition metal. The present study of 
graphene may give a key to understanding the proper- 
ties of other carbon nanostructures, such as carbon nan- 
otubes. Up to now, atomic scale simulations are inca- 
pable of explaining the chiral-selective growth observed 
experiment ally j 37 i 38 The approach proposed here could 



help to identify individual healing mechanisms during 
growth that produces perfect tube structures and those 
favoring a definite chirality. 
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